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ABSTRACT 

The  atmospheric  correction  of  thermal  infrared  (TIR)  imagery  involves  the  combined  tasks  of  separation  of  atmospheric 
transmittance,  downwelling  flux  and  upwelling  radiance  from  the  surface  material  spectral  emissivity  and  temperature. 
The  problem  is  ill  posed  and  is  thus  hampered  by  spectral  ambiguity  among  several  possible  feasible  combinations  of 
atmospheric  temperature,  constituent  profiles,  and  surface  material  emissivities  and  temperatures.  For  many  materials, 
their  reflectance  spectra  in  the  Vis-SWIR  provide  a  means  of  identification  or  at  least  classification  into  generic  material 
types,  vegetation,  soil,  etc.  If  Vis-SWIR  data  can  be  registered  to  TIR  data  or  collected  simultaneously  as  in  sensors  like 
the  MASTER  sensor,  then  the  additional  information  on  material  type  can  be  utilized  to  help  lower  the  ambiguities  in  the 
TIR  data.  If  the  Vis-SWIR  and  TIR  are  collected  simultaneously  the  water  column  amounts  obtained  form  the 
atmospheric  correction  of  the  Vis-SWIR  can  also  be  utilized  in  reducing  the  ambiguity  in  the  atmospheric  quantities. 
The  TIR  atmospheric  correction  involves  expansions  in  atmospheric  and  material  emissivity  basis  sets.  The  method  can 
be  applied  to  hyperspectral  and  ultraspectral  data,  however  it  is  particularly  useful  for  multispectral  TIR,  where  spectral 
smoothness  techniques  cannot  be  readily  applied.  The  algorithm  is  described,  and  the  approach  applied  to  a  MASTER 
sensor  data  set. 
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1.  INTRODUCTION 

Space-based  Thermal  IR  (TIR)  imagery  is  a  powerful  remote  sensing  method  that  can  operate  in  nighttime,  fog  and  haze 
conditions.  We  include  in  the  thermal  infrared  all  emission  dominated  radiance  in  the  MWIR  and  LWIR  between  4-15 
pm.  Several  approaches  for  enhancing  the  information  content  of  TIR  imagery  have  been  studied  in  recent  years, 
including  increasing  the  spectral  resolution  to  the  hyperspectral  or  ultraspectral  level.  High  spectral  resolution  is  valuable 
for  retrieving  atmospheric  properties  and  for  separating  surface  temperature  and  emissivity.  Currently,  ambiguities  in 
separating  temperature,  emissivity  and  atmospheric  effects  complicate  the  interpretation  of  the  imagery  using  multiple 
wavelength  bands.  High  spectral  resolution  TIR  requires  a  larger  spatial  footprint  for  an  equivalent  signal-to-noise  than 
multispectral  TIR.  The  large  footprints  make  the  detection  of  small  objects  more  difficult.  The  approach  of  fusing 
multispectral  TIR  imagery  with  short-wave  data  helps  overcome  these  problems  and  has  several  other  important 
practical  advantages.  Vis-SWIR  data  can  be  collected  ahead  of  time,  even  several  days  in  advance,  and  from  a  different 
platform  than  the  TIR  sensor’s.  The  Vis-SWIR  data  must  be  collected  during  the  day,  and  such  sensors  generally  have 
very  good  signal-to-noise  and  spatial  resolution.  TIR  provides  for  both  day  and  night  observations. 

The  material  reflectance  spectra  in  the  Vis-SWIR  provide  a  means  of  identification  or  at  least  classification  into  generic 
material  types,  vegetation,  soil,  etc.  If  Vis-SWIR  data  can  be  registered  to  TIR  data  or  collected  simultaneously,  as  with 
sensors  such  as  the  MASTER  sensor,  then  the  additional  information  on  material  type  can  be  utilized  to  help  lower  the 
ambiguities  in  the  TIR  data.  In  this  effort  we  apply  FLAASH1  to  Vis-SWIR  to  obtain  reflectance  maps.  These  are  used 
to  identify  materials  by  matching  scene  reflectances  with  library  reflectances. 

To  ensure  that  the  scene  reflectances  correspond  more  closely  to  single  materials,  we  use  SMACC,2  an  end-member 
code,  to  determine  the  unique  pixel  spectra  in  the  scene.  Identification  of  a  material  type  is  not  always  sufficient  to 
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determine  the  reflective  properties,  particularly  of  natural  materials.  Various  states  of  health  and  leaf  cover  lead  to  a 
wide  variety  of  vegetation  reflectance  spectra.  Vegetation  spectra  in  the  Vis-SWIR  are  used  to  remotely  determine 
vegetation  conditions.3  These  variations  extend  into  the  TIR.  The  use  of  basis  sets  which  include  these  variations  have 
been  applied  to  successfully  unmix  pixel  spectra.5’6  If  the  Vis-SWIR  and  TIR  are  collected  simultaneously,  the  water 
column  amounts  obtained  from  the  atmospheric  correction  of  the  Vis-SWIR  can  also  be  utilized  to  reduce  the  ambiguity 
of  the  atmospheric  quantities  in  TIR  imagery. 

In  this  paper  we  address  multispectral  TIR  atmospheric  correction,  a  critical  step  in  the  fusion  process.  Several  efforts 
have  been  made  to  perform  atmospheric  correction  on  hyperspectral  TIR  data.  These  include  the  in-scene  atmospheric 
correction  method,  ISAC,7which  finds  scaled  estimates  of  atmospheric  transmittance  and  path  radiance.  It  must  be 
applied  to  blackbody  pixels,  water  or  high  emissivity  vegetation  areas.  The  method  requires  a  statistical  sampling  of  high 
emissivity  pixels  over  a  range  of  ground  temperatures,  making  it  difficult  to  apply  to  all  scenes.  For  low  emissivity 
materials  in  the  corrected  scene,  a  source  term  is  derived  that  includes  contributions  from  reflected  downwelling 
atmospheric  radiance. 

Several  methods  use  forward  modeling  including  a  look-up  table  (LUT)  approach  based  on  a  grid  of  MODTRAN 
calculations8  and  spectral  smoothness9.  The  approach  provides  a  complete  separation  of  the  scene  into  pixel  emissivities, 
temperatures  and  atmospheric  profiles.  Material  emissivities  are  spectrally  smooth  compared  to  gas  phase  atmospheric 
absorption  and  emission  features.  Broad  channels  in  multispectral  sensors  limit  the  utility  of  the  smoothness  criterion  for 
such  sensors.  Other  forward  models  use  grids  of  calculations  to  form  bases  utilizing  component  analysis.  ITernandez- 
Baquero  and  Schott10  have  developed  an  atmospheric  correction  model  based  on  canonical  correlation  analysis,  CCA. 
This  method  correlates  a  basis  of  atmospheric  transmittance,  path  radiance  and  downwelling  radiance  to  the  scene 
radiances  to  select  an  atmosphere,  and  is  followed  by  a  temperature  emissivity  separation  step.  The  method  shows 
promise,  however  the  correlation  coefficients  coupling  the  three  atmospheric  quantities  to  total  radiance  are  different  and 
the  atmosphere  selected  does  not  necessarily  correspond  to  a  single  profile.  Villeneuve  and  Stocker11  have  developed 
OPRA,  an  oblique  projection  technique  to  perform  atmospheric  correction.  Using  a  grid  of  atmospheric  transmittances 
and  path  radiances  computed  from  MODTRAN,  OPRA  uses  two  oblique  subspace  projections  to  estimate  the 
atmospheric  path  radiance  and  transmittance  separately.  Independent  projections  can  lead  to  the  transmittance  and  path 
radiance  being  modeled  in  different  regions  of  the  subspaces  and,  hence,  the  predicted  atmospheric  quantities  need  not 
correspond  to  a  single  atmospheric  profile.  It  is  currently  implemented  for  blackbody  or  highly  emissive  materials. 
Flealey  and  Slater12  use  library  reflectance  spectra  and  generate  a  grid  of  radiances  for  a  material  for  a  wide  range  of 
atmospheric  conditions  and  illumination.  Basis  sets  formed  from  these  radiances  span  subspaces  for  atmosphere  and 
environment  independent  material  detection. 

The  basis  of  our  approach  is  a  forward  model  that  involves  creation  of  a  set  of  radiances  that  span  the  space  of  the  scene 
pixel  radiances  and  use  a  physically  constrained  least  squares  approach  to  simultaneously  estimate  the  atmospheric 
quantities  and  the  ground  material  emissivity  and  temperature.  Atmospheric  downwelling  effects  are  included.  The 
estimated  atmospheric  quantities  are  traceable  to  a  single  atmospheric  profile.  The  methodology  is  described  in  Section 
2.  Initial  applications  of  the  method  are  described  in  Section  3. 

2.  TIR  ATMOSPHERIC  CORRECTION  ALGORITHM 

The  radiance  model  is  briefly  described  and  then  the  atmospheric  correction  strategy  is  developed.  Throughout  this  text 
we  suppress  the  spectral  index  for  simplicity.  Flowever,  each  of  the  radiative  transport  terms  is  dependent  on  the  spectral 
wavelength  interval.  The  quantities  that  we  generate  are  the  mean  emissivity  and  the  mean  emission  of  a  pixel.  We 
consider  the  most  general  case  of  a  mixture  of  materials  in  a  pixel  at  a  variety  of  temperatures.  In  this  case  the  mean 
emission  is  not  due  to  a  single  emissivity  blackbody  product,  but  rather  from  a  sum  of  such  products. 

The  model  used  for  the  radiance  at  the  sensor,  L  ,  of  a  mixed  pixel  is 

L  =  Erm  +  (l  —  e)D  +  U 
and  the  radiance  at  the  source,  S,  is  given  by 


(1) 


S  =  (L-U)/x 


(2) 


where  D  is  the  downwelling  radiance  that  is  transported  to  the  observer,  U  is  the  upwelling  or  path  radiance,  T  is  the 

direct  ground  to  observer  transmittance,  e  is  the  mean  pixel  emissivity  and  Em  is  the  mean  pixel  source  emission.  The 
mean  emissivity  of  a  pixel  containing  r  materials  is  given  by 


e  =  ll£Pap 

P 


(3) 


where  S  is  the  spectral  emissivity  of  material  p  and  up  is  the  abundance  of  material  p  in  the  pixel.  The  abundances 
r 

sum  to  unity,  X  a p  —  1  ■  It  follows  that  the  mean  source  emission  is  given  by 
P 

Em  =  YJspapYJfp,B(Tl) 

p  i  (4) 

where  /p  i  is  the  fraction  of  material  p  that  is  at  temperature  Th  The  temperature  distributions  of  each  material  must  also 
sum  to  unity,  X  fp  /  =  /  •  B(T)  is  the  Planck  blackbody  function.  There  are  a  total  of  r  materials  at  5  temperatures  in 
the  FOV. 

This  model  includes  all  of  the  atmospheric  interactions  of  importance  for  the  TIR.  The  downwelling  radiance  is 
composed  of  indirect  scattered  radiance  and  atmospheric  thermal  emission.  The  model  currently  ignores  reflected 
emissions  from  other  ground  sources,  a  component  that  will  be  added  to  the  model  in  the  future.  The  total  hemispherical 
reflectance,  p  =  1  -e  ,  is  used  in  place  of  the  full  BRDF  (Bi-directional  Reflectance  Distribution  Function).  Bi-directional 
reflectance  effects  are  small  in  the  TIR  spectral  region  where  emission  dominates.  The  effects  of  directional  emissivity 
can  be  included  into  the  pixel  source  emission  term.  Directional  effects  will  be  more  important  when  processing  data  sets 
collected  at  off  nadir  angles. 

The  basic  problem  of  atmospheric  correction  in  the  TIR  is  the  uncoupling  of  the  atmospheric  and  material  contributions, 
that  is,  atmosphere  profde,  emissivity  and  temperature  separation.  Strong  coupling  among  the  quantities,  together  with 
the  sensitivity  of  the  ground  temperature  predictions  and  lower  altitude  atmospheric  temperature  predictions  to  small 
variations  in  emissivity,  make  extraction  of  the  of  the  atmospheric  and  the  ground  quantities  very  difficult.  Initial 
guesses  to  any  of  these  quantities  strongly  influence  the  predicted  results  for  the  other  two  quantities.  We  developed  a 
solution  via  a  triple  expansion,  with  simultaneous  extraction  of  atmospheric  quantities,  emissivities  and  temperatures. 

A  radiance  basis  is  constructed  by  combining  atmospheric  and  material  emissivity  and  temperature  basis  sets.  The 
atmospheric  basis  is  obtained  by  running  SMACC  on  a  look-up  table  formed  from  a  set  of  MODTRAN  calculations  that 
tabulate  the  transmittance,  upwelling  and  downwelling  radiance,  that  we  refer  to  as  the  TUD,  for  a  grid  of  atmospheric 
representations  {tk,  Uk,  Dk}.  The  index  k  designates  position  within  the  spanned  atmospheric  parameter  space.  A  basis 
set  of  material  emissivities,  {em}  can  be  selected  from  library  or  previously  processed  data.  The  emissivity  basis  should 
be  selected  to  span  the  anticipated  ground  emissivities  {sp}.  Finally,  a  set  of  temperatures,  {Tn}  can  be  selected  that  span 
the  anticipated  ground  temperatures,  and  a  basis  of  blackbody  emissions  is  formed  using  the  Planck  radiance  expression. 
A  radiance  basis  is  constructed  from  the  three  sets  by  combining  all  possible  combinations  of  atmospheric  profile, 
emissivity  and  temperature.  The  set  of  radiances,  { },  is  given  by 


^  k,m,n  -  emB(Tn  Jlk  +  (1  em  )^k  +  ^ k 


(5) 


Formally,  the  pixel  radiance  can  be  estimated  by  calculating 


min  ||  L  £  A k,m,n Qc  ,m,n 

k,m,n 


(6) 


with  respect  to  the  coefficients  C/(mn ,  subject  to  physical  constraints,  i.e.  by  solving  a  constrained  least  squares 

problem.  The  radiance  basis  is  constructed  from  the  TUD  basis,  the  emissivity  basis  and  a  grid  of  temperatures  using 
Equation  5  rather  than  performing  MODTRAN  runs  for  each  profile,  emissivity  and  temperature.  This  leads  to 
significant  reductions  in  processing  time,  since  each  MODTRAN  calculation  needs  only  to  be  performed  once  per 
profile,  and  not  repeated  for  each  emissivity  and  temperature  combination. 

Once  the  expansion  coefficients  have  been  found,  the  estimates  of  the  atmospheric  transmittance,  downwelling  and 
upwelling  path  radiance  are  given  by 

t=  ,m,n  >  D  =  Y  DkCk,m,n  an^  U  =  Y^! kC k ,m,n  ■  respectively.  (7) 

k,m,n  k,m,n  k,m,n 


The  expansion  coefficients  are  the  same  for  all  three  quantities,  allowing  the  three  atmospheric  quantities  to  be  derived 
from  a  single  profile.  If  the  coefficients  are  constrained  to  be  positive,  the  resulting  combination  of  the  basis  provides  an 
interpolation  of  the  basis  atmospheres  as  the  estimate  to  the  unknown  atmospheric  TUD.  Other  estimated  quantities 
include  the  emissivity,  source  emission,  apparent  blackbody  emission  and  apparent  temperature. 

The  emissivity  of  the  pixel  can  be  obtained  from  the  reflected  downwelling  contribution. 


S  = 


k,m,n 


,m,n 


S  k.m,n 

k,m,n 


(8) 


It  is  a  weighted  average  of  the  basis  of  emissivities  used  in  the  calculations.  The  ground  material  emission  is  determined 
from 


Em  = 


Y  em  B(Tn  Jlk  Qc  ,m,n 
k,m,n _ 

Ytk^k  ,m,n 
k,m,n 


(9) 


If  it  is  further  assumed  that  either  a  single  material  is  present  within  the  footprint  of  the  pixel  or  that  all  materials  within 
the  footprint  are  at  the  same  temperature,  the  emission  can  be  factored  into  an  apparent  Planck  blackbody  function. 


B 


app  ~ 


Y  em  B(Tn  X k  Q:  ,m,n 
k,m,n _ 

YemBSk,m,n 

k,m,n 


(10) 


Also,  an  additional  expression  for  the  ground  emissivity  can  be  obtained  as 


(11) 


Z  em  B(  Tn  )tk  Q:  ,m,n 
_  k,m,n _ 

Z  B(  Tn  )t/c  Cfc mn 

k,m,n 

The  two  emissivity  estimates,  Equations  8  and  11,  should  agree  closely  for  pixel  ground  materials  at  the  same 
temperature.  The  apparent  temperature  is  obtained  as  a  function  of  the  brightness  temperature  of  the  apparent  blackbody 
emission. 

Refined  estimates  to  the  atmospheric  and  ground  quantities  can  be  made  based  on  a  finer  grid  of  basis  functions  centered 
around  the  current  output  estimates,  or  on  additional  criteria  such  as  emissivity  spectral  smoothness,  or  smoothing  the 
apparent  blackbody  estimate  by  selecting  a  single  temperature  from  the  spectrally  dependent  brightness  temperatures 
retrieved  form  the  apparent  blackbody  function  ,  Equation  10. 

Implementation  of  the  simultaneous  approach  leads  to  a  large  basis  set.  For  example,  if  there  are  25  atmospheric  basis 
functions,  4  emissivity  basis  functions  and  2  temperatures,  then  the  total  number  of  radiance  basis  functions  is  25x4x2 
=200.  The  basis  will  generally  be  redundant  and,  as  written,  the  formal  solution  is  underdetermined.  The  rank  or  the 
number  of  linearly  independent  basis  spectra  will  be  less  than  the  number  of  radiance  basis  spectra.  This  means  that 
there  are  many  possible  solution  vectors.  The  most  that  can  be  simultaneously  taken  as  independent  equals  the  number 
of  spectral  bands,  but  the  actual  number  may  be  far  fewer  since  the  spectra  are  highly  correlated  for  these  wavelengths. 

There  are  several  approaches  that  can  be  taken  to  solve  the  underdetermined  constrained  equations.  We  seek  an 
approach  that  will  allow  us  to  obtain  estimates  of  the  atmospheric  and  material  properties.  Here  we  consider  the 
combined  use  of  constraints  and  a  singular  value  decomposition  (SVD)  basis.  This  leads  to  a  new  set  of  equations  that 
is  not  underdetermined.  We  form  a  matrix, A  of  radiances  whose  columns  are  the  {A kmn)-  Singular  value 

decomposition  can  be  applied  to  the  A  matrix.  The  expansion  can  be  performed  using  the  significant  singular  vectors, 
those  with  large  singular  values,  rather  than  the  A  basis  directly.  The  singular  vectors  with  non-zero  singular  value  form 
a  linearly  independent  set.  The  SVD  is  given  as 


A  =  WsVT  .  (12) 

The  columns  of  the  matrix  W  are  the  left  singular  vectors,  the  basis  spectra  that  will  be  used  in  the  expansion.  The 
matrix  s  is  square  and  diagonal  and  contains  the  non-zero  singular  values  on  its  diagonal,  and  the  matrix  V  contains  as  its 
columns  the  right  singular  vectors.  The  new  problem  is  to  find  the  vector  x  whose  elements  satisfy 

min  1 1 L  -  X  WfXf  \  |  (13) 

i 

subject  to  the  constraints  on  the  elements  of  x  that  contributions  of  the  redundant  A  basis  remain  positive.  The 
constraints  placed  on  x  can  be  expressed  in  terms  of  the  back-transformation 

Vs~]x  =  C  .  (14) 

The  constraint  relations  are 

Vs~JX  >  0  .  (15) 

To  solve  the  underdetermined  problem  posed  by  Equation  6,  we  solve  the  fully  determined  Equation  13,  then  back- 
transform  via  Equation  14  to  obtain  the  set  of  coefficients,  Cfc  m  n  . 

The  individual  coefficients,  Cjcm  n  are  not  well  estimable  due  to  the  high  correlation  among  the  columns;  however,  the 
linear  combinations  of  the  columns  that  lie  in  the  subspace  of  the  dominant  singular  vectors  are  estimable.13 


The  set  of  coefficients,  C/(  m  n  derived  here  are  unique  in  that  they  are  all  positive  or  zero  and  they  are  the  components 
of  the  solution  vector  of  shortest  length.14  The  SVD  approach  with  constraints  leads  to  a  physically  possible  solution. 
We  use  the  back-transformed  coefficients  to  estimate  the  unknown  atmospheric  and  material  quantities  as  outlined 
above. 

An  alternative  approach  to  determining  a  reduced  basis  is  to  use  an  end-member  procedure  to  find  a  linearly  independent 
set  of  vectors  in  the  radiance  basis.  In  the  work  described  here  we  used  a  hybrid  approach,  using  SMACC  to  find  a  basis 
to  represent  an  atmospheric  look-up  table,  then  performing  a  SVD  on  the  radiance  basis  constructed  with  the  reduced 
atmospheric  basis,  emissivity  basis  and  a  grid  of  temperatures. 

3.  APPLICATION  OF  THE  TIR  ATMOSPHERIC  CORRECTION  METHOD 

The  initial  testing  of  the  algorithm  included  synthetic  radiances  constructed  from  a  loam  emissivity  at  a  specified  ground 
temperature  and  a  selected  atmosphere  from  a  LUT.  The  selected  atmosphere  was  not  among  those  selected  as  a  basis  to 
represent  the  LUT,  rather  it  was  the  atmosphere  that  was  least  well  represented  by  the  selected  basis.  The  synthetic 
unknown  is  used  to  evaluate  the  ability  of  the  algorithm  to  extract  atmospheric  and  material  quantities  for  situations  in 
which  the  basis  of  atmospheres,  emissivities  and  temperatures  span  or  approximately  span  the  simulated  radiance  and  for 
cases  where  one  of  the  quantities  is  not  well  represented.  Given  the  constraints,  a  radiance  spectrum  will  be  in  the  span 
of  our  basis  if  the  atmospheric  quantities  are  represented  by  a  positive  linear  combination  of  TUD  basis  members,  the 
emissivity  is  represented  by  a  weighted  average  of  the  emissivity  basis  and  the  ground  temperature  lies  within  the  range 
of  selected  temperatures. 

An  atmospheric  transmittance,  upwelling  and  downwelling  radiance  (TUD)  look-up  table  was  generated  for  the 
MASTER  sensor15  mid-wave  and  long-wave  IR  channels.  The  atmospheric  profile  characteristics,  water  vapor  column 
amounts  and  profile  temperatures,  and  sensor  altitude  were  selected  so  that  the  TUD  would  be  appropriate  for  processing 
a  MASTER  California  agriculture  data  cube  that  had  been  processed  in  the  Vis-SWIR  with  FLAASH. 

The  atmospheric  temperature  profile  model  requires  four  temperatures  and  altitudes  to  specify  it,  see  Figure  1.  The  top 
two  altitudes  and  temperatures  are  above  the  sensor  altitude  of  3.18  kilometers.  The  altitudes  for  T |,  T2,  T3  and  T4  were 
selected  as  Okm,  2km,  7km  and  13km,  respectively.  An  altitude  range  can  be  specified  for  which  the  temperature  profile 
remains  constant  at  T3,  but  in  this  case  it  was  specified  at  a  point.  All  profiles  in  the  TUD  had  the  same  upper 
atmospheric  profile  with  temperature  values  T3=255K  and  T4=215K.  An  equally  spaced  grid  of  atmospheric 
temperatures  at  the  ground,  T ,,  ranging  from  280K  to  310K,  was  selected.  The  equally  spaced  second  temperature  grid 
T2  at  2  kilometers  was  selected  to  range  from  268K  to  288K.  The  grid  of  water  vapor  column  amounts  for  the  TUD 
ranged  from  70%  to  100%  of  a  standard  humidity  profile  that  is  similar  to  the  mid-latitude  summer  profile  contained  in 
MODTRAN.  A  basis  to  represent  the  TUD  was  obtained  by  using  SMACC  to  determine  the  endmembers  of  the  TUD. 
The  endmember  basis  is  determined  sequentially  with  an  improved  representation  of  the  TUD  at  each  step.  At  each  step 
the  atmosphere  which  is  least  well  approximated  by  the  basis  is  added  to  the  basis.  A  total  of  twenty  five  endmembers 
were  found  to  adequately  represent  the  remaining  TUD  members. 

The  synthetic  unknown  radiance  spectrum  was  constructed  from  an  atmosphere  in  the  TUD.  The  selected  atmosphere 
was  not  part  of  the  basis.  It  was  selected  as  the  one  which  was  least  well  approximated  by  the  basis.  This  selection 
should  test  the  ability  of  the  basis  to  model  any  atmosphere  within  the  convex  span  of  the  LUT’s  atmospheres.  The  lower 
atmosphere  temperature,  Tb  of  this  atmosphere  was  the  hottest  in  the  LUT,  Ti=310K.  The  value  of  T2=276K  was  a 
mid-value  in  the  LUT.  The  water  vapor  content  was  the  lowest  possible  in  the  LUT  table,  a  value  of  70%  of  the  standard 
profile.  This  particular  atmosphere  was  excluded  from  consideration  for  the  basis.  A  loam  spectrum  from  the  ASTER 
library16  was  selected  to  represent  the  ground  material.  Its  library  number  is  87P4453c  and  it  is  described  as  a  “very  dark 
grayish-brown  loam”.  This  spectrum  was  excluded  from  the  emissivity  basis  sets,  although  five  other  loam  spectra  were 
used  to  form  a  loam  emissivity  basis.  The  ground  temperature  for  the  synthetic  unknown  radiance  spectrum  was  set  at 
300K,  a  full  10K  below  the  lower  atmospheric  temperature,  Tb 


Atmosphere  Profile 


Temperature 


Figure  1.  Definition  of  temperature  profile  for  atmospheric  TUD  calculations. 

The  radiance  basis  for  the  calculations  was  constructed  from  the  endmember  basis  of  the  TUD  together  with  various 
combinations  of  emissivity  bases  and  ground  temperature  grids.  As  a  preliminary  test,  we  used  the  same  material 
emissivity  in  the  basis  as  in  the  unknown.  This  let  us  determine  the  adequacy  of  the  atmospheric  endmembers  to  describe 
the  unknown  atmosphere  and  the  feasibility  of  using  a  grid  of  ground  temperatures  to  estimate  the  true  ground 
temperature.  We  also  performed  calculations  using  a  spanning  basis  of  other  loam  emissivities  that  did  not  contain  the 
actual  ground  emissivity  spectrum  and  a  non-spanning  basis  that  consisted  of  vegetation  spectra.  The  results  showed  the 
sensitivity  of  the  atmosphere  and  temperature  determination  to  emissivity  estimates.  The  results  for  the  non-spanning 
emissivity  basis  lead  to  unphysical  atmospheric  quantities  and  poor  ground  temperature  estimates.  The  results  for 
incorrect  guesses  of  ground  temperature  lead  to  poor  estimates  of  emissivity  and  poor  or  unphysical  estimates  of  the 
atmospheric  quantities.  The  results  for  a  grid  of  two  temperatures,  which  bracketed  the  actual  ground  temperature,  were 
as  good  as  guessing  the  temperature  exactly,  and  the  ground  temperature  could  be  extracted  from  the  brightness 
temperature  of  the  estimated  blackbody  function. 

The  results  for  several  of  the  initial  calculations  are  summarized  in  Figure  2.  Calculations  labeled  “truth”  in  Figure  2 
correspond  to  the  inputs  to  the  test  radiance.  The  atmospheric  basis  was  used  for  all  other  calculations.  Calculations 
labeled  “Full  Model”  use  the  emissivity  basis  that  includes  other  loam  spectra  and  a  pair  of  ground  temperatures,  298K 

and  302K  which  bracket  the  test  ground  temperature  of  300K.  Calculations  labeled  “Model  8  known”  use  the  same 
emissivity  spectrum  as  the  test  calculations  and  the  pair  of  ground  temperatures,  298K  and  302K  which  bracket  the  test 
ground  temperature  of  300K.  Calculations  labeled  “Model  T=298K”  and  “Model  T=302K”  use  the  basis  of  ground 

emissivities  and  single  ground  temperatures  which  are  too  low  and  too  high,  respectively.  Calculations  labeled  “Model  £ 
=  vegetation”  used  an  emissivity  spectrum  of  vegetation  and  the  pair  of  ground  temperatures  of  298K  and  302K.  The 

closeness  of  the  “Model  £  known”  results  to  the  “Truth”  indicate  that  the  atmospheric  basis  does  span  the  atmospheres 
in  the  LUT  and  that  adequate  ground  temperature  estimates  can  be  made  by  selecting  a  sufficient  grid  of  ground 
temperatures  as  model  input.  The  agreement  of  the  “Full  Model”  with  the  “Truth”  in  Figure  2  gives  an  indication  of  the 
agreement  that  can  be  expected  for  favorable  conditions  of  a  spanning  atmospheric  basis,  a  spanning  ground  emissivity 
basis  and  a  pair  of  bracketing  temperatures.  The  remaining  calculations  illustrate  the  effects  of  selecting  a  model  ground 
temperature  too  low  or  too  high,  or  of  selecting  the  ground  emissivity  of  the  wrong  material  type.  If  the  model 
temperature  is  selected  too  low,  298K,  an  incorrect  atmosphere  with  a  higher  low  atmospheric  temperature  Tj  results, 
and  the  predicted  emissivity  is  too  high  and  has  an  incorrect  shape.  If  the  model  temperature  is  too  high,  302K,  an 
incorrect  atmosphere  with  a  lower  low  atmospheric  temperature,Ti  results  and  the  predicted  emissivity  is  too  low  and 
has  the  wrong  shape.  If  an  incorrect  emissivity  spectrum  is  selected,  vegetation  in  place  of  loam,  the  overall  shapes  of 
the  predicted  atmospheric  quantities  are  in  disagreement  with  the  true  quantities  and  the  ground  temperature  is  estimated 
as  299. 5K. 
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Figure  2.  Results  of  calculations  estimating  the  emissivity  (a),  transmittance  (b),  upwelling  radiance  (c)  and  downwelling  flux 
(d)  compared  with  the  synthetic  unknown  radiance  inputs  (Truth). 


Application  to  the  MASTER  Scene 

The  first  application  we  made  of  the  multispectral  TIR  atmospheric  correction  approach  to  real  data  was  to  a  MASTER 
sensor  California  agricultural  scene  illustrated  in  Figure  3.  FLAASH  was  used  to  determine  reflectances  in  the  Vis- 
SWIR  as  well  as  the  humidity  and  visibility.  SMACC  was  used  to  determine  endmember  reflectances  in  the  VIS-SWIR 
and  the  reflectance  endmembers  were  matched  against  library  spectra  to  identify  material  type.  The  endmember  material 
abundances  were  used  to  determine  groups  of  spectra  of  the  same  type.  The  groups  with  large  numbers  of  pixel  spectra 
were  chosen  for  atmospheric  correction.  Nine  endmember  groups  were  selected.  These  groups  were  primarily  varieties 
of  vegetation  with  some  soil  with  one  group  being  primarily  soil. 

The  emissivity  basis  used  for  the  calculations  included  four  spectra:  two  green  vegetation  spectra,  one  of  high  emissivity 
and  one  of  low  emissivity,  a  dry  grass  spectrum  and  a  soil  spectrum.  The  vegetation  emissivity  spectra  were  taken  from 
those  made  available  by  Wan17.  The  spectra  illustrate  the  variability  of  reflectance  properties  of  vegetation  in  the  TIR 
spectral  region.  The  soil  spectrum  is  the  average  spectrum  of  25  ASTER  library  soil  emissivity  spectra.  The  spectra  are 
illustrated  in  Figure  4. 


Figure  3.  MASTER  California  agricultural  scene. 

An  initial  range  of  bracketing  ground  temperatures  was  chosen  between  290K  and  310K,  then  narrowed  to  296K  and 
304K.  The  ability  to  successfully  use  the  linear  approximation  to  unmix  TIR  data  has  been  recently  illustrated.18 
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Figure  4.  Emissivity  Basis  for  calculations  include  a  high  emissivity  chlorophyll  containing  vegetation  spectrum,  a  low  emissivity 
chlorophyll  containing  vegetation  spectrum  ,  a  dry  grass  spectrum  and  a  soil  spectrum. 


Two  atmospheric  basis  sets  were  utilized,  the  one  used  in  the  application  to  the  synthetic  spectrum,  and  a  second  smaller 
basis  derived  from  an  atmospheric  LUT  in  which  the  humidity  was  restricted  to  3085.5  atm-cm,  the  value  retrieved  from 
the  Vis-SWIR  spectral  region  using  FLAASH.  In  this  smaller  LUT  the  upper  atmospheric  temperatures  were  held  fixed 
and  the  lower  temperatures  were  varied  from  290K  to  310K  for  Tj  and  280K  to  300K  for  T2.  Differences  among  the 
resulting  atmospheric  corrections  obtained  from  the  two  atmospheric  basis  sets  illustrate  the  advantages  of  fully  utilizing 
results  of  atmospheric  correction  of  Vis-SWIR  data  to  aid  in  the  TIR  atmospheric  correction. 

Results  for  the  two  sets  of  calculations  applied  to  the  nine  endmember  groups  are  reported  in  Figures  5  and  6  for  the 
atmospheric  basis  with  variation  in  atmospheric  humidity  and  for  the  atmospheric  basis  with  humidity  restricted, 
respectively.  Differences  between  the  two  sets  of  calculations  include  a  much  narrower  spread  in  ground  temperatures 
but  a  broader  spread  in  lower  atmospheric  temperature  for  the  variable  humidity  atmospheric  basis.  Both  humidity  and 
lower  atmospheric  temperature  variations  contribute  to  the  variability  among  predicted  atmospheric  transmittance,  and 
upwelling  and  downwelling  radiance.  The  ground  temperature  and  the  downwelling  radiance  estimations  are  the  most 
sensitive  to  the  small  differences  in  extracted  emissivity.  These  are  the  quantities  that  are  the  most  different  for  the  two 
sets  of  calculations.  For  high  emissivity  materials  the  downwelling  contribution  to  the  radiance  is  small,  and  moderate  to 
large  variations  in  the  downwelling  radiances  lead  to  very  small  changes  in  total  radiance.  The  use  of  the  auxiliary 
information  obtained  from  FLAASH  about  the  atmosphere  leads  to  a  greater  relative  consistency  for  the  predicted 
atmospheric  quantities  among  the  set  of  endmember  pixels. 
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Figure  5.  Results  for  atmospheric  correction  of  nine  pixel  groups  using  the  variable  humidity  atmospheric  basis,  (a)  transmittance, 
(b)  upwelling  radiance  (c)  downwelling  flux,  (d)  emissivity  (e)  ground  temperature. 
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Figure  6.  Results  for  atmospheric  correction  of  nine  pixel  groups  using  the  fixed  humidity  atmospheric  basis,  fixed  to 
match  the  prediction  of  FLAASH,  (a)  transmittance,  (b)  upwelling  radiance  (c)  downwelling  flux,  (d)  emissivity  (e)  ground 
temperature. 

CONCLUSIONS 

Our  approach  to  TIR  atmospheric  correction  uses  a  forward  model  with  bases  of  atmospheric  TUDs,  emissivities  and 
blackbody  emissions.  It  is  a  triple  expansion  with  assumptions:  (1)  the  true  atmospheric  transmittance,  upwelling  and 
downwelling  radiances  are  positive  linear  combinations  of  the  bases  of  the  {tk,  Uk,  Dk}  derived  from  MODTRAN,  and 
(2)  all  quantities  related  to  the  ground  material(s)  are  weighted  averages  of  the  equivalent  basis  set  quantities. 
Mathematically,  these  assumptions  are  equivalent  to  assuming  that  the  atmospheric  quantities  are  in  the  convex  cones  of 
the  corresponding  atmospheric  basis  quantities  and  that  the  ground  material  quantities  are  in  the  convex  hull  of  the 
derived  ground  material  bases.  To  derive  a  pixel  temperature,  the  additional  assumption  that  the  pixel  contains  only  one 
material  or  the  materials  in  the  pixel  are  at  the  same  temperature  is  required. 

The  process  of  atmospheric  correction  and  emissivity/temperature  separation  of  TIR  data  is  ill-posed.  The  use  of  a 
minimum  residual  criterion,  common  to  this  method  and  other  subspace  projection  techniques,  is  hampered  by  shallow 
minima,  with  a  variety  of  combinations  of  compensating  ground  and  atmospheric  conditions  leading  to  close  fits  to  the 
data.  The  atmospheric  and  ground  material  quantities  are  strongly  coupled  and  estimating  them  is  sensitive  to  small 
variations  in  emissivity  and  atmospheric  profile.  Auxiliary  information  and  constraints  make  important  contributions  to 
reduce  the  ambiguities.  Knowledge  about  the  ground  materials  allows  for  smaller,  but  more  relevant  emissivity  basis 
sets.  The  Vis-SWIR  data  can  be  used  to  identify  material  types  and  to  select  pixel  groups  for  TIR  atmospheric 
correction.  The  variability  of  the  reflectance/emissivity  of  natural  materials  in  the  Vis-SWIR  and  the  TIR  spectral 
regions  can  be  treated  through  the  use  of  bases  of  reflectances/emissivities.  Additional  information  about  atmospheric 
conditions  can  further  reduce  the  parameter  space  required  for  the  atmosphere  description. 

The  method  can  be  refined  by  repeating  calculations  using  finer  basis  grids  after  initial  estimates  are  made.  We  are 
subjecting  the  method  to  further  testing  against  synthetic  data  and  we  plan  to  apply  it  to  TIR  data  for  which  ground  truth 
is  available. 
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